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A clear look at fuzzy PI control 



By Peter Morgan 



The second 
installment of 
InTech's advanced 
control tutorial 
series looks at one 
use of fuzzy logic in 
process control. 



Fuzzy logic for indusrrial automation is gain- 
ing popularity, but how can it be applied to con- 
tinuous process control? One way is with fuzzy PI 
control, a technique for adjusting the output 
from a proportional-plus-integral (PI) algorithm 
based on the results of an associated fuzzy logic 
algorithm. 

Before considering fuzzy PI in detail, it is 
instructive to review the conventional PI algo- 
rithm in its incremental (or rate) form. This form 
of the algorithm is commonly used where the 
controller interfaces to a servo motor for final ele- 
ment positioning. For purposes of this evalua- 
tion, it allows the characteristics of fuzzy PI to be 
more easily compared with conventional PI. 

The proportional-plus-integral algorithm in 
its series form is: 
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= manipulated variable (controller output) 
= controller or proportional gain 
= integral action time (sec) 
= controller error = +/- (process variable - 
set point), with + for direct-acting loops 
and - for reverse action 
Presuming final output manipulation is an 
integration either in the program or at the final 
element — for example, when the output is treat- 
ed as velocity demand for a servo motor — it is 
necessary to differentiate the PI algorithm to 
retain Pi-type adjustment of the manipulated 
variable. In other words, 
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If the control system samples at an interval of 
T seconds, then 
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where 

Amv = required change in manipulated 

variable at end of each sample interval 



Ae = change in error between samples 
Solving for the change in manipulated variable, 
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If, after every scan, the output changes by 
Amv as defined, then the output adjustment will 
match that obtained from a conventional PI algo- 
rithm implemented in a continuous (analog) 
controller. 

Position demand can be obtained by integrat- 
ing Amv through a simple accumulation: 
mv = mv lasI + Amv 

Since the conventional, or crisp, control algo- 
rithm determines the change in output from the 
error and change in error, the fuzzy algorithm 
should rely on input memberships for these same 
variables. Consider three categories, or grades, for 
the error and Aerror: large-positive (L PO s)> OK, 
and large-negative (L NEG ). As few as two member- 
ship grades can specify the input variables, but three 
grades more clearly demonstrate the principle. 

Ignoring for the moment how the input vari- 
ables are to be treated quantitatively, we can 
establish a rule set from our understanding of the 
incremental form of the PI algorithm. The first 
rule might be: 

IF {E = L POS } AND (AE = L POS ) THEN 
{AMV=L POS } 

where E and AE are the scaled error and change 
in error, and AMV is the scaled change in output. 

This is entirely consistent with the equation 
for the rate form of the PI controller. It is also 
consistent with a rule based on intuition. In other 
words, if the error is already large and is increas- 
ing, then act to increase the output a lot. 

The second rule could be: 
IF {E = L POS } AND (AE = L NEG } THEN 
(AMV = OK} 

In other words, if the error is large and posi- 
tive, but changing to be less positive, then do 
nothing (much!). This too might have been 
derived from the rate form of the conventional 
algorithm in so far as e and Ae would tend to 
cancel, resulting in little change in the output. 
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hand, the process loop has two comparable first- 
order lags in the control valve and the transmitter. 
Also, on average, the process-loop dead time using 
control valve C is significantly longer when revers- 
ing direction. With the latter type of process 
response, it is unlikely that high controller gain 
can be used in a conventional proportional-plus- 
integral control scheme. 

10% steps 

Figure 6 indicates that with significandy larg- 
er amplitude steps, the motion of control valve A 
was velocity limited; it followed a nearly straight 
line, indicating the actuator being filled with the 
positioner relay wide open. The response of con- 
trol valve B was still asymmetrical, but noticeably 
faster. The response of control valve C was faster 
than in the 2% steps, but still had an exponential 
behavior. Due to its single-stage positioner, con- 
trol valve C was still significantly slower than the 
other designs. 

As the step size in the open-loop tests becomes 
progressively larger, the control engineer should 
consider whether the control system is actually 
capable of delivering such a large command 
change to the final control element. 

100% steps 

The 1 00% steps were conducted in a separate 
test with the pump operating and the same load 
valve position as in the previous tests. The sample 
rate in the latter test was 55 per second. Full- 
range tests are sometimes called stroking-time 
tests. In most applications, this data is relevant 
only for emergency conditions. 

Figure 7 shows that, for these largest steps pos- 
sible, control valve A became significantly slower 
than the others. The major factors causing the 
slow stroking are the larger volume and lower 
supply pressure of this style of diaphragm actua- 
tor. The relays in the positioner of control valve A 
also had smaller flow area than that of control 
valve B. If fast response is desired for very large 
steps, use of a volume booster between the posi- 
tioner and actuator would be appropriate with 
control valve A. The booster action should be 
limited to large-amplitude input changes, using a 
needle valve bypass around the booster relay. 

The response of all three control valves was 
velocity limited over part of the 100% stroke. 
Control valve B still had faster response than con- 
trol valve C. This apparendy was due to a larger 
limiting airflow area in the positioner. 

It is essential to recognize the strong amplitude 
dependence of the response-time data. The user 
must distinguish between the small-amplitude 
and large-amplitude requirements of the process 
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Figure 7. Stem-position response for 100% steps ("stroking-time" test) 



control loop. Stroking-time data (Figure 7) can- 
not be scaled down to predict the small-amplitude 
behavior (Figure 5) of the control valve. 

Try this at home 

Users can easily apply step-response tests to 
obtain their own open-loop control valve perfor- 
mance data. They can measure dead band by 
making small moves about an operating point, 
and determine dynamic response. However, 
always use a step size appropriate to the control 
objective. Using steps larger than the distur- 
bances the control loop will experience may 
result in selection of inappropriate equipment. 
Measurement of process variable changes, under 
process loading, is the ultimate confirmation of 
the desired behavior. 

It is clear that broad performance claims, 
without considering all the components of the 
control valve, are not supported by the data. 
Superior performance requires superior design 
and integration of the valve, actuator, and posi- 
tioner. For example, a belief that piston actuators 
always give better performance is clearly erro- 
neous; combining a piston actuator with a low- 
performance positioner creates a slow-moving 
control valve. On the other hand, a low-friction 
valve coupled with a conventional diaphragm 
actuator and a two-stage positioner can give 
excellent dead band and small-amplitude 
response time. IT 
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The complete rule set for the PI algorithm is 
shown in Figure 1, where for convenience the 
rules are shown on a grid (or matrix). 



Error 




Figure 1. Fuzzy PI controller will change manipulated variable 
based on a "rule matrix." 

Normalizing £ and Ae to bring E and AE 
within a range from -1 to +1 enables transforma- 
tion of actual values to membership grades. The 
-1 to +1 range is used in most discussions of fuzzy 
sets, but, in practice, an alternative range may 
better suit the particular control system. For 
example, an integer value of 0-4,095 is more 
appropriate for a PLC that does not support 
floating-point math. 

E = ^ 



Assume raw error is in percent of transmitter 
range (a benefit when comparing scaling factor 
values with conventional gain and integral action 
settings). Then, requiring maximum incremental 
adjustment of the output when the raw error is 
20% and change in raw error is only 1 % per scan 
leads to e max =20 and Ae max =l. Since e max and 
Ae max determine the degree of action for a given 
error and change in error, these scaling factors 
provide a means of tuning the algorithm. It will 
be shown later how these factors relate quantita- 
tively to settings in the conventional algorithm. 

As E, or AE, varies between -1 and +1 mem- 
bership functions, such as those in Figure 2, 
determine the degree with which it qualifies in 
each of the membership categories. 

In practice, the grades can be easily computed 
using simple algebra: 
Grade(OK) = l-abs(E) 
Grade(L NEG ) = abs(E) for E less than or equal 
to zero 

= for E greater than zero 
Grade(L POS ) = abs(E) for E greater than or 
equal to zero 
= for E less than zero 

For illustrative purposes, assume control error 
is 5% of transmitter range and change in error 
per scan is 0.2%. If the error and Aerror maxi- 
mum values are 30% and 0.5%, then E and AE 
are 0.167 and 0.4, respectively. The grades, or 




Grade value in OK 
category 

Grade value in L P0C 
category 



degree of membership, of E and AE in the OK 
category are thus 0.833 and 0.6, respectively. 
Figure 3 shows how other grades follow from the 
application of the membership function equa- 
tions for LpQ S and L NE q 



Figure 2. Fuzzy membership functions for 




Figure 3. Membership grades for 5% control error and 0.2% 
change in control error per scan 



From the rule matrix for fuzzy PI (Figure 1), 
when E is Lp s and AE is L PO s the change in 
output is required to be a large increase (L INC ); 
that E and AE are only partially L PO s means that 
the change in output for this rule must be only 
partially "large increase." 

For Boolean variables, defined by the states 1 
or 0, it is advantageous for manipulative purpos- 
es to consider the "AND" function as the multi- 
plication of two, or more, logical states. In other 
words, 1x0 = 0, lxl = 1, 0x0 = 0. Taking the 
minimum Boolean state will also yield the 
desired logical result. So, too, in fuzzy logic is the 
minimum of the two input grades equivalent to 
the resultant grade for the consequent output. 

Considering the grades for E 



and AE in the L POS category 
(Figure 3) and the PI rule set 
(Figure 1), we would determine 
the output grade as 0.167 Lj NC . 
Applying this principle to each PI 
rule results in output member- 
ships shown in Figure 4. 

Referring to the rule matrix 
(Figure 1), we may say, for in- 
stance, that the output will mod- 
erately increase if any of the antecedent rules call 
for a moderate increase; this is equivalent to an 
"OR" function. Noting that taking the maxi- 
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Figure 4. Moderate changes, M (NC and 
M DEC , each have two grades, and no 
change has three grades in the output 
membership matrix. 
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Figure 5. Obtain crisp output by finding 
the center of area of fuzzy output mem- 
bership function. 
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Figure 6. Representing output member- 
ship functions as series of discrete val- 
ues simplifies calculation of crisp output 



mum of Boolean states provides 
us with the correct logical result 
for all combinations of inputs for 
the OR function, for a given out- 
put membership category we can 
take the maximum of the evalu- 
ated grades. Figure 5 shows the 
tabulated grades and resultant 
output membership functions for this input 
condition. 

Determining the incremental, or crisp, output 
requires determining the center of area (or center 
of gravity) of the combined output membership 
functions. This is most easily 
accomplished by representing 
each membership function by a 
number of vertical elements. 
Figure 6 shows a single member- 
ship function represented by nine 
elements, including zero-value 
elements. Element value will 
equal the minimum of the value 
determined by the functions tri- 
angular envelope or the particular category's 
grade. 

For each output category, element values may 
be mapped onto the output axis, as Figure 7 
shows. In practice, this is accomplished by calcu- 
lating V;, the value of the function at the zth posi- 
tion on the output axis from: 
V^Max^Jv.Fj) 

where Fj through Ft are the membership func- 





tion values applying at the rth position. 

In this case, ignoring zero grades, only two 
membership functions apply at any position. At 
position +3, for instance, OK (value 0.25) and 
M lNC (value 0.4) apply; there- 
fore, at position +3 the output 
function value V ; is 0.4. 

Referring to Figure 7, if V; is 
the value at the kh interval on 
the Amv axis, then the center of 
area can be evaluated by dividing 
the sum of the individual 



moments about the zero axis by the sum of the 
individual values multiplied by the coefficient of 
normalization, 8 for this case: 

12 

IVixi 

AMV = ^ 

-12 

For this example, the center of area, AMV, will be 
0.327. 

The actual change in output is the per unit 
change in output multiplied by the output scal- 
ing factor, or: 
Amv = AMV x Amv max 

Thus, if the scaling factor, Amv max , is 2.5, 
when the algorithm output is 1, the change in 
final element position demand will be 2.5%. 
While it is possible to choose a scaling factor that 
matches the maximum rate of change in output 
to the velocity limit of the output actuator, this 
has an ameliorating effect on the response. For 
most purposes, it is better to use Amv max =100 
so that proportional-type action is unrestricted. 

To obtain the final position demand, the cal- 
culated change in output is summed with the last 
output: 

mv = mv^ + Amv 

For this example, with an output scaling fac- 
tor of 100, the output change would be 32.7%. 

The complete process, described schematically 
in Figure 8, continuously repeats this single-cycle 
computation (performed once per controller 
scan) with outputs calculated using the error and 
Aerror obtained at the beginning of each scan. 

For study and prototyping purposes, a spread- 
sheet can easily implement this computational 
approach. Representing the axis of Figure 7 with 
a series of 25 spreadsheet cells, each cell will 
incorporate the formula for the membership 
function value V;. 

Pi/fuzzy PI action comparisons 

The first step in comparing the conventional 
and fuzzy PI characteristics is evaluating the 



Figure 7. Discrete approximation to out- 
put membership functions map to output 
"axis" in process of determining crisp 



Figure 8. Fuzzy proportional plus integral 
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incremental output for the complete input range 
of inputs. Plotting incremental output on a verti- 
cal axis against error and Aerror on orthogonal 
horizontal axes yields a three-dimensional surface 
representing incremental output. Surface plots 
for both algorithms (Figure 9) show that, while 
the surfaces are similar, the fuzzy algorithm sur- 
face is characterized by cusping that results from 
the nonlinear relationship between change in 
output and error or Aerror. The ratio of fuzzy 
outputs to conventional outputs for each error 
condition, plotted in Figure 10, more clearly 
reveals the difference between actions. 




anal algorithr 



Figure 10. Plotting fuzzy output/co 
more clearly illustrates difference between actions over input 
range. 

From Figure 10, we may deduce that the fuzzy 
algorithm introduces nonlinearities that result in 
nonlinear dependencies in both proportional and 
integral actions. For small inputs, the fuzzy algo- 
rithm applies more aggressive correction than 
conventional PI. When Aerror has polarity oppo- 
site that of the error, i.e., when the process is 
moving in the right direction, the algorithm 
tends to be less aggressive. 

Pi/fuzzy PI parameter equivalence 

While effective gain and integral action time 
for fuzzy PI control exhibit nonlinear depen- 



dence on error and change in error, it is possible 
to establish the approximate equivalence of fuzzy 
PI parameter settings to the gain and integral 
action settings for the conventional algorithm. 

For the fuzzy algorithm, when error is zero, a 
change in error equal to the maximum value 
adopted in the scaling (normalization) process 
gives rise to a 50% change in output (for Amv max 
= 100). 

More generally: 

AmV max 

Amv = 0.5 x 



Ae , 



■x Ae 



Comparing this with the proportional com- 
ponent of the conventional PI algorithm in its 
incremental form, it follows that: 

AmV max 



G = 0.5x 



Ae , 



Also, when the change in error is zero, an 
error equal to the maximum value adopted in the 
scaling process gives rise to a 50% change in out- 
put (for Amv max =100). In other words, 

AmV max 

Amv = 0.5 x x e 

£ max 

Referring to the integral term in the incre- 
mental form of the conventional algorithm, we 
may deduce that: 

T AmV max 

G X — = 0.5 X 

T £ max 

Rearranging and substituting for G, the 
equivalent integral action time can be obtained 
as: 

T = T x - ° 



Ae , 



Tuning the fuzzy algorithm to meet the 
adopted response criteria for the process model 
used in this study requires setting Ae max = 65 and 
£ max = 2,500. These settings transform to an 



Terminology 



mv = 



e = 



x = 



G = 



Ae= 



T = 



c max> 
Ae max= 



AMV = 



controller output, 
i.e., manipulated 
variable 

controller error 

integral time = 1/1, 
where I is the inte- 
gral action rate 

controller or 
proportional 



Amv m , y = 



i in error 
during sampling 

sampling interval 



error and Aerror, 
which cause maxi- 
mum incremental 
adjustment of 



maximum incre- 
mental adjustment 
of controller output 

scaled error = 



scaled incremental 
change in manipu- 
lated variable = 

Amv 



Amv. 
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Figure 11. 
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Fuzzy PI and conventional PI responses for 10% load disturbance are indistinguishable. 




Figure 12. Fuzzy PI exhibits slightly less 
overshoot than conventional PI for 10% 
set-point change with controller opti- 
mized for load disturbance. 



equivalent gain of 0.8 and integral action time of 
19.2 seconds for a sample time of 0.5 seconds. 
The conventional PI controller requires gain of 
1.25 and integral action time of 19.25 seconds. 
In other words, to obtain similar response, inte- 
gral action times closely correspond, but the gain 
for conventional PI is somewhat higher than the 
transformed value. This is consistent with the 
demonstrated nonlinearity of the fuzzy algorithm 
and the fact that controllers are tuned to meet 
response objectives for a 10% load disturbance. 

Who's the better performer? 

In comparing responses, we strived to estab- 
lish a good reference response for conventional 
PI. Often, tuning the reference PI goes no further 
than applying approximate methods, such as 
Ziegler-Nichols. In other words, without fine- 
tuning the reference algorithms against specified 
criteria, definitive response comparisons cannot 
be made. 

The process model used in this analysis com- 
prises a dead time approximation based on three 
3-second lags and a 20-second first-order lag; 



actuator response is incorporated as a 5%/sec rate 
limit. 

We tuned the conventional PI algorithm so 
that, for a load disturbance, integral error squared 
is minimized without there being any overshoot 
in the simulated process variable, pv. In other 
words, response is critically damped. We tuned 
the fuzzy algorithm againsr the same objective. 
The resulting responses (Figure 11), are barely 
distinguishable from one another. 

Since load disturbance response is the same, 
we might ask if there is a difference in response 
for set-point changes. 

For set-point changes (Figure 12), the fuzzy 
algorithm exhibits slightly lower overshoot and 
less tendency for overshoot to increase as set- 
point changes become smaller. However, com- 
paring performance indexes, evaluated as inte- 
gral-error squared, shows the conventional algo- 
rithm does somewhat better. We may quite fairly 
conclude that, for Pi-type algorithms, there is no 
substantive difference in performance between 
the conventional algorithm and fuzzy algorithm. 

Then why fuzzy logic? 

While fuzzy logic may be no better than the 
conventional PI algorithm, it has comparable 
performance and can have its place. 

Fuzzy logic lends itself to evolution of algo- 
rithms based on observations of required control 
actions typically carried out by human operators. 
In particular, the equivalence of fuzzy operations 
to those used in Boolean logic greatly assists in 
transforming operating practices or observations 
into deterministic control actions. 

Where processes are typically controlled by 
the operator with some success, careful evalua- 
tion of operator methods can lead to the estab- 
lishing of rule sets that can be used directly in the 
formulation of a fuzzy algorithm, which can be 
applied to control the process, reduce the require- 
ment for constant operator intervention, and 
bring about performance improvements through 
the consistent application of the best operating 
practice. IT 
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